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Braneworld models in extra dimensions can be tested in laboratory by the coupling of the radion 
to the Standard Model fields. The identification of the radion as a canonically normalized field 
involves a careful General Relativity treatment: if a bulk scalar is responsible for the stabilization 
of the system, its fiuctuations are entangled with the perturbations of the metric and they also have 
to be taken into account (similarly to the well-developed theory of scalar metric perturbations in 4D 
cosmology with a scalar field). Extracting a proper dynamical variable in a warped geometry /scalar 
setting is a nontrivial task, performed so far only in the limit of negligible backreaction of the scalar 
. . , field on the background geometry. We perform the general calculation, diagonalizing the action up 

' to second order in the perturbations and identifying the physical eigenmodes of the system for any 

' amplitude of the bulk scalar. This computation allows us to derive a very simple expression for the 

, exact coupling of the eigenmodes to the Standard Model fields on the brane, valid for an arbitrary 

(N , background configuration. As an application, we discuss the Goldberger-Wise mechanism for the 

' stabilization of the radion in the Randall-Sundrum type models. The existing studies, limited to 

, small amplitude of the bulk scalar field, are characterized by a radion mass which is significantly 

. below the physical scale at the observable brane. We extend them beyond the small backreaction 

\^ ' regime. For intermediate amplitudes, the radion mass approaches the electroweak scale, while its 

\ coupling to the observable brane remains nearly constant. At very high amplitudes, the radion 

mass instead decreases, while the coupling sharply increases. Severe experimental constraints are 
T-H ' expected in this regime. 

> : 
o\ . 

OO ■ I. INTRODUCTION 

' Braneworlds embedded in higher dimensions have opened new interesting directions in high energy particle physics, 

T^lj- , general relativity, and cosmology. The Randall-Sundrum (RS) model with warped 5D bulk geometry between two 

■ plane-parallel orbifold branes jlj has been particularly studied. Its low energy spectrum contains a scalar field, denoted 
as radion Q, which is associated to the inter-brane modulus. The radion is coupled to Standard Model (SM) fields 

• living on one brane, which gives very interesting phenomenology for accelerator experiments ^J. The coupling is of 

\^ ^' the form 

^: a„* = |t;, (1) 

J> I where is the trace of the energy momentum of the SM fields, Q refers to the radion field, while the coupling A is 
' close to the electroweak scale Jl, |3 0, • 

. The radion originates from the scalar fluctuations of the bulk geometry. In the RS model the two branes are at 
' indifferent equilibrium: the interbrane distance is a free parameter, corresponding to a massless radion. The situation 
is changed if a bulk scalar field is introduced to stabilize the system, as in the Goldberger-Wise mechanism 
The identification and the study of the properties of the radion field are less straightforward in this case. The initial 
studies, including the original work 0, were based on the effective 4D theory, obtained integrating an approximate 
ansatz for the geometry and the scalar field profile along the extra coordinate. However, the introduction of the 
bulk scalar makes the gravity /scalar system more complex, and requires a self-consistent treatment based on the 5D 
Einstein eqs. In particular, the scalar fluctuations $ of the bulk metric are sourced by the fluctuations Scf) of the 

bulk scalar field, and the coupled linearized system (<&, (50) also has to be treated self-consistently. The 5 dimensional 
equations for the system were derived in the laasic papers 0, ||| . The analysis of |^ is particularly advanced and clear 
with respect to the coupling of the modes, so we will often refer to it for comparison. The properties of the radion 
obtained in these studies are nonetheless in good agreement with the ones from the 4D effective theory (the agreement 
is actually much better than what has been claimed sometimes in the literature, as we will discuss below). 

All the existing studies share one limitation: it is always assumed that the bulk scalar field has a negligible 
backreaction on the background geometry (the only effect being that it lifts up the flat direction corresponding to the 
radion). Phenomenological studies have so far been restricted to this case, first because it is only in this limit that 
the bulk geometry is (approximately) AdS as in the RS proposal, but also because it is the only case for which the 
properties of the radion were known. Indeed, in this limit the contribution of 54> to the wave function of the radion 
is negligible, and the radion behaves approximately like the metric perturbation $ 'Sj. As a consequence, the exact 
diagonalization of the coupled ($, 5(j}) system in not required in this limit. 
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In contrast to these works, the aim of the present paper is to study the scalar perturbations for the general situation 
in which the bulk scalar </) has arbitrary amplitude. The main motivation is that very small </) should not be considered 
as completely natural, since the idea of RS is to avoid strong hierarchies besides the one given by the warp factor. ^ 
There is no incompatibility between a sizable tf) and the main idea of : although the bulk geometry will be in general 
different from pure AdS, one can still naturally achieve a strong warping at the observable brane, as needed to solve 
the hierarchy problem. In particular, this is also true for the Goldberger-Wise mechanism ;3] and its generalizations. 
The main obstacle is that the mixing between the different scalar perturbations cannot any longer be neglected, and 
the coupled system of perturbations has to be diagonalized. This study is still missing in the existing literature, and it 
is the main result of the present work. As we will show, the coupled ($, 5(j)) pair contains a single dynamical variable 
which we denote by v, and which is a linear superposition of the pair. The equation for v reduces to an oscillator-like 
equation. More importantly, the action of the system can be shown to be diagonal in terms of the 4D Kaluza-Klcin 
modes of v , which we denote by Qn . In this way, we obtain the exact 4D action for the physical scalar modes of 
the theory. The action differs from the effective (approximate) 4D theory which has been computed in the literature 
by neglecting the Stj) modes, and - moreover in general - by starting from a 5D ansatz which does not satisfy the 
linearized equations for the perturbations (the effective action is derived and discussed in Appendix ICl to emphasize 
the difference from the exact 4D theory discussed in the rest of the paper) . 

The final result can be easily summarized. Decomposing the metric perturbation in terms of the 4D KK modes, 
$ = ^^j^ $„ (y) Qn (x) {y denoting the extra coordinate), we find that the quadratic action for the scalar modes can 
be cast in the diagonal form 

S = Y.Sn = Y.Cn f d^xQn[n-ml]Qn , (2) 

n n 

where are the physical masses obtained from an eigenvalue boundary problem in the extra dimension (see below). 
As we remarked, the radion corresponds to the mode with the lowest mass. The fields Qn are scalars fields of the 
exact 4D theory, and their decoupled actions can be trivially quantized. The coefficients C„ are given by 



_ 3 m3 ry« dy 







3M3 



(3) 



In this expression, A is the warp factor, M the fundamental scale of gravity in 5 dimensions, and prime denotes 
differentiation with respect to y . This result is particularly simple and completely general, since it holds for arbitrary 
background configurations A and cj) . The derivation, to be discussed in detail below, exploits the strong similarity of 
the system with the one of scalar metric fluctuations in the scalar field inflationary cosmology. This analogy allows 
us to identify the dynamical variable v, which in inflationary cosmology is known as the Mukhanov-Sasaki variable 
[id lll| , and which is the starting point of our computation. 

The requirement that the modes Qn are canonically normalized, C„ = 1/2 , fixes the normalization of , which in 
turn defines the strength of the coupling of the radion and of the other KK modes to the SM fields on the observable 
brane (the computation and diagonalization of the action is indeed necessary for this task. The normalization of the 
modes cannot be derived from the Einstein equations for the perturbations, which are all linear in the perturbations). 
Similarly to the initial expression (Q, the coupling A„ of the radion and the other KK modes to brane fields is given 
by 



^ _ Qn rpfj. 

•t-int.n — . 



1 




(yo) 


A„ 




2 



(4) 



where yo is the location of the observable brane. Due to the explicit knowledge of ©, we can determine the coupling 
of all the perturbations for arbitrary background configuration. Our general calculation confirms the accuracy of the 
results of ^6J in the limit of small (p . ^ Our results extend these calculations to a general amplitude of the field </> and 
to an arbitrary geometry A . 

The paper is organized as follows. In section^lwe present the action and the relevant equations of the system. The 
main result summarized above is derived in section lTTll while it is applied to a couple of specific examples in section llVl 
one of which - first proposed in Q - is a generalization of the Goldberger-Wise mechanism 3] for the stabilization of 



^ A larger is also more natural (at least for the Goldberger-Wise mechanism) since it increases the fraction of initial conditions for which 

the stabilization is achieved, if the radion is initially displaced from the stabilized value. See Appendix |C] for a discussion. 
^ We thank C. Csaki for discussions and clarifications on this issue. 
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the radion of the Randall- Sundrum model [ij. The existing studies of the phenomenology of this model, limited to 
small amplitude of the bulk scalar field, are characterized by a radion mass which is significantly below the physical 
scale at the observable brane. This hierarchy decreases as (j> increases, and, fore sizable values of cj) , the radion 
mass approaches the electroweak scale. On the contrary, the coupling of the radion to the observable sector remains 
constant for a wide range of . The radion mass reaches a maximal value for a finite value of . At higher amplitudes 
the mass decreases, while its coupling to the SM sharply increases. In section we verify numerically some of these 
analytical results: we follow the oscillations of a two-brane system around a static configuration with the numerical 
code developed in 12]; the frequency spectrum of the oscillations exhibits resonances precisely in correspondence to 
the masses of the modes given by the analytical computation. The higher resonances are excited by the their coupling 
to the radion in the nonlinear regime. Once their amplitude becomes small, they become constant in time, confirming 
the fact that they are decoupled at the linear level. The results of the paper are discussed in section IVll where we 
also comment on more general related issues. For example, we discuss preheating effects due to the oscillations of the 
radion, and the linearized 4D gravity on the branes. The paper is concluded by three Appendices. In the first one 
we compute the quadratic action for the perturbations. Appendix IbI presents analytical results in the small limit. 
Finally, Appendix ICl contains the computation of the effective 4D action for the radion, and the comparison with the 
exact results discussed in the rest of the paper. 



II. SET-UP 



We start from the action 

rvo 



S ■ 







E 

'i=l,2 



(5) 



where the first term is the bulk action, and the sum contains contributions from the two branes, located at y = and 
y — y^) = — 2/0 along the extra dimension. The overall factor of 2 multiplying the first term accounts for the integration 
over the intervals [— 2/o,0] and [0,?/o], which are identified by an S'^/Z2 symmetry. To shorten the notation, the bulk 
cosmological constant is included in the bulk potential V , while brane tensions are included in the brane potentials 
U . The quantity [K] denotes the jump of the extrinsic curvature at the two branes, while M is the fundamental scale 
of gravity in 5 dimensions. We choose the following line element 



ds^ = A [yf [(1 + 2^-) ?7^^ dx^" dx" + (1 + 2$) dy^] 



(6) 



where Greek indices span the 1 + 3 coordinates parallel to the flat branes, and rj — diag (— 1, 1, 1, 1) . Namely, we 
have conformal coordinates for the warped background geometry, while the scalar metric perturbations 5* and $ are 
expressed in the (generalized) longitudinal gauge. ^ Perturbations of the bulk scalar field are denoted ip {x, y) = 
4>{y) + 54>{x,y) . We are interested in static warped configurations (and, in addition, flat branes), so that the two 
background quantities A and 4> depend on the extra coordinate y only. Concerning the background quantities, 
variation of the action |(SJ) gives the bulk equations (here and in the following, prime denotes derivative wrt y ) 



0" + 3^0' 
^ - 0— 



^2 y = , 

^'2 



3M3 



.A_ 
^2 



0'2 



A^V 



as well as the junction conditions 



l2 



U 

6M3 



A 



±- 



U' 



(7) 



(8) 



We need to consider a smaller set of scalar metric perturb atio ns with respect to the one discussed in brane cosmology, when the induced 
metric on the brane is of FRW type, see for instance Il3l . When, as in our case, the brane is maximally symmetric, the se other 
modes are decoupled from the ones we consider here, and they are actually components of the massive gravitational waves |14| . Scalar 
perturbations also include the disp lac ements £,i of the two branes along the y coordinate. However, in absence of fluids or kinetic terms 
for the scalar on the brane, = Il5l . so that we do not need to include them in the calculation. In Section llVl we consider brane fields, 
to compute their interaction with the radion and the KK modes. However, it is accurate to treat them as probe fields (consistently with 
the interaction picture), which do not modify the setting discussed here. 



4 



where the upper and lower sign refers to the brane at ?; = and y = i/q , respectively (prime on the potentials V and 
U denotes differentiation with respect to 4>). 

As with the background, not all of the equations for the perturbations are dynamical, reflecting the fact that not 
all the modes given above are physical. The system of linearized Einstein equations for the perturbations include the 
two constraint equations 



A, 



3M3 A' 



= 



as well as a dynamical equation. By using the constraint equations, the latter can be written as 



□ - 



jl/z)" 
dy^ l/z 



U = 



where we have defined 



A' 



U 



^3/2 



(9) 



(10) 



(11) 



and where □ is the D'Alambertian operator in 4 d. 
at the two branes, 

(A^ 



Finally, there are two junction conditions for the perturbations 



3M3 ' 

±-AU"i 
2 



(12) 



The first of these conditions could have also be obtained from the two constraint equations (O evaluated at the two 
branes. The second condition instead contains additional, physically relevant information, relating 5(f> to the "mass 
parameters" of the scalar field at the two branes. As can be intuitively expected, very high U" force Scj) ^ {) at the 
two branes. This is the case which is typically considered in the literature (starting from the original proposal 3]), 
since the results simplify significantly in this limit. In the next section we perform the identification and quantization 
of the perturbations for arbitrary [/", and we give general formulae valid for arbitrary bulk and brane potentials for 
(j) . The examples of section IIVI are instead restricted to high U" . 



III. QUANTIZATION OF THE ACTION FOR PERTURBATIONS AND THE IDENTIFICATION OF 

THE PHYSICAL MODES. 

From the analysis of the previous section, we see that scalar metric perturbations are coupled to the fluctuations of 
the bulk scalar field responsible for the stabilization of the system. The covariant formulation |(SJ fully underlines the 
symmetries of the theory. However - as in all gauge theories - the price to pay is the introduction of auxiliary degrees 
of freedom, so that it is not straightforward to identify and decouple the physical modes of the system. However, this 
is very important if we want to know their properties, namely their masses and their coupling to the (electroweak) 
sector on the observable brane. The mass spectrum for the perturbations can be obtained directly from the equations 
of motion. The easiest way to do so is to first identify some linear combination of the perturbations for which the 
dynamical equation can be cast in a Schrodinger-type form 0, |^ 0, . This is the case for the variable U defined 
in eq. (|llll . although other combinations can be used (see below). After solving this equation, the spectrum of the 
system can be obtained from the boundary conditions (|12|l at the two branes. 

This procedure, however, does not identify which combination v of the perturbations corresponds to the physical 
degrees of freedom, nor does it determine the proper normalization of the modes, since the bulk and boundary 
conditions arc linear in the perturbations. ^ On the other hand, this information is needed to determine the coupling 
of the physical modes to the observable sector, see eq. To obtain it, one has to diagonalize the action 52 for 



* On a technical level, W is a combination li = cJJi + C2W2 of the two linearly independent solutions of equation IIUI . The boundary 
conditions at the two branes allow the determination of the ratio ci /c2 , as well as the physical mass of the mode (defined as the 
eigenvalue of □ , see below) , but they do not allow to determine ci and C2 independently. 



5 



these modes, which is computed by expanding the initial action ^ up to second order in the perturbations. This is a 
more difficult exercise than simply obtaining the equations of motion for the perturbations, and indeed it is the main 
original result of the present work. 

We divide the discussion in two subsections. In the first one we quantize the system for the case in which the bulk 
scalar field (j) is absent. This is a particularly simplified situation, since in this case there is only one physical scalar 
perturbation. This computation, which is the only one given in the literature, allows us to determine the normalization 
of the radion when the bulk scalar field is absent or when its backreaction on the bulk geometry can be neglected ^ . 
As shown in JJj, the contribution of Scj) to the wave function of the radion is negligible in this limit, and the coupling 
of the radion to brane fields is (approximately) given by the coupling of $ . In subsection llll BI we perform the general 
calculation, valid for a bulk scalar field of any amplitude. This computation confirms the accuracy of the analysis 
of in the limit in which the backreaction of (j) is small and can be treated perturbatively. In addition, it gives new 
interesting results. First, as we mentioned, it provides the coupling of the radion to brane fields for arbitrary values 
of . Second, it allows the computation of the coupling of the other KK modes, which cannot be approximated by 
the analysis of llll Al even in the small </> limit. The reason for this is the discontinuity in the number of modes with or 
without the bulk scalar field: when (p is absent, there is only one physical scalar degree of freedom (the radion). On 
the contrary, there is a whole tower of KK modes when the scalar field is present (no matter how small its amplitude 
is). Hence, the study of the action without (f> does not provide any information about the KK modes. 



A. Computation without the bulk scalar field 



It is very instructive to study what happens when the scalar field is absent. Expanding the action (jSJ at quadratic 
order in the perturbations, and imposing ^E* — — $/2 , we find 



S2 [f^O]-- 
From eq. Q we have 



2Ar 



12 



.A^ 
A^ 



A' 

$' + 2— $ = ^ $cxA"2 
A 



(13) 



(14) 



Hence, in this case there is only one scalar eigenmode, corresponding to the radion field, while the tower of KK modes 
is absent 0| . We can then eliminate <&' in favor of <& , so that the action (|13|) becomes 



52 0] = 2M^ / d^xA^ 



A/2 Al. 

^^,y-''a,$a.$ + 3(2— --)$ 



A 



(15) 



The second term vanishes due to the background equations {T)). We thus recover the well known result of a massless 
radion in the absence of any bulk field. Decomposing 



$(x,2/) = $(2/)Q(x) , 



the above action becomes 



3Af^ 



dyA^^' \ rfW9^Qa,Q = -- / d''xr,^"'^f,Q^^.Q 



(16) 



(17) 



The last condition (i.e. the request that the 4D field Q is canonically normalized) specifies the normalization of $ , 
which in turns determines the coupling of the radion with the fields on the observable brane. The quantization of the 
last expression is straightforward. 

When the bulk scalar is absent, the function Q , which encodes the physical degree of freedom, obeys the free wave 
equation DQ = . The free bulk wave equation is satisfied by the bulk gravitational waves Hab- The 5 bulk degrees 
of freedom of Hab are projected to the 4D brane as two tensor modes (corresponding to usual 4D gravitons), two 
gravi-vectors (which are however absent at orbifold branes) and one gravi-scalar. When the bulk scalar is absent, we 
suggest interpreting the massless radion Q as the massless gravi-scalar projection of the bulk graviton. 



^ The action is obtained by setting (f) = Sep = in the general expression I18i given below. 
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When the branes are stabihzed by the bulk scalar field, the massless gravi-scalar projection is absent However, 
this degree of freedom does not disappear. The dynamical degree of freedom will be in this case a linear combination 
of the metric perturbation $ and the fluctuation of the scalar field 50 , and it emerges as a massive radion at the 
brane, rather than as a gravi-scalar (see the next section for details). The situation is similar to the re-shuffling of the 
degrees of freedom which occurs during a symmetry breaking, where here the symmetry breaking is associated with 
the stabilization of the intcrbrane distance. 



B. General computation with a bulk scalar field 



The expansion of the action @ at second order in the perturbations gives 



S2 = 2Ar / (Tx { J] 



A' 



dfiS(f)diyS(l) 



+ *) 9,* 



1 



A'' 



+121- VP' + 2VI;) _ I ^ - 121^ j f -$2 + ) - 



d^^A^ 



2L/*2 -I- 2C/'(5(/)* + -V'Sfj)^ 



2R'P 
1 



i'2 + 6,I,'2 



— ^',50' ($ - 4*) - 



(18) 



The derivation is given in Appendix Due to the presence of the bulk scalar field and its fluctuations, the general 
expression (jl8|l is more involved than and the quantization is more difficult. There are however similarities with 
the well known problem of scalar perturbations in 4D inflationary cosmology with a scalar fleld, which fortunately can 
be exploited in the case at hand. ^ In the cosmological set-up, the inflaton field is the analogous of the bulk scalar (j) , 
and the time t is analogous to the bulk coordinate y (the scale factor a (t) is the analogue of the warp factor A{y)). 
The evolution equations for the infiationary case constitute an initial value problem (initial conditions specified at 
some early time during infiation), which is in our case replaced by a boundary value problem (junction conditions at 
the two branes). The analogy cannot be pursued up to the quantization (in both cases, a canonical quantization is 
performed, and the t ^ y analogy has to be abandoned). However, it allows us to identify the dynamical variable v 
which is the starting point for the quantization of (|18|) . There are two complications compared to the infiationary case. 
The first is that, from the AD point of view, the system (|18|l contains an infinite KK tower of scalar perturbations, 
while there is only one scalar mode in the infiationary case (as it was also the case in the previous subsection). The 
second is the presence of two branes as boundaries of the y coordinate. Total derivatives in y cannot be dropped 
(unlike the cosmological situation, where the space has no boundaries). As we shall see, they are actually crucial to 
ensure the orthogonality of the different modes in the action. 

Starting from the action H18|) . we first eliminate $ by use of the first of lO (nondynamical equations, have to be used 
for the quantization. They compensate for the presence of non-physical degrees of freedom in the set of perturbations). 
We then identify total derivatives, and simplify the final expression by use of the background equations. Finally, we 
use '3/ = — $/2 to express (and its derivatives) back in terms of $ . Although the calculation is rather involved, the 
final result is particularly simple 



d^xdy 



□ 



dy^ 
3 A^ 
A A' 



z 
z 



A' 



$ 2 



~A 



A^V ^A U"(j)' 



The variable z was defined in Hll() . The combination 

$ A' 



V = z 



A 



3M 



43 



{zuy 



(19) 



(20) 



is the 5D generalization of what is commonly denoted as Mukhanov-Sasaki variable [Tollll|. and it is the dynamical 
variable in the bulk. It is gauge invariant, i.e. it does not change under infinitesimal redefinitions of the coordinate 



The reduction to the dynamical action that we give here extends the one performed in llSlligl for the case of cosmological perturbations. 
Although we show it here in longitudinal gauge, we have actually derived it in arbitrary gauge for the perturbations. 
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system. The bulk equation for v is immediately obtained from the action H19|l . and it can be shown to be equivalent 
to the dynamical eq. H1U|) expressed in terms of U 



□ 



w = 



(21) 



As we mentioned, the boundary terms cannot be neglected. Their presence enforces the hermiticity of the Lagrangian 
operator. Starting from only the bulk term, we get '' 

<f z 



1 



(fx { Vl 



□ 



(fx 



V2 - V2 



dy2 



Vl 



, , A' A^V AU"<i)' 



(22) 



The last term precisely cancels with the difference of the boundary terms in (|19|) . This guarantees the orthogonality 
of the KK modes of the scalar perturbations. The different modes arise from the decomposition 



F(a:,y) =^F„(y) Q„ (x) , 



(23) 



where F is any of Stf), $, and v (the same functions Qn have to be used in order to satisfy the coupled equations H12|) 
and 1(201) • The equations for the perturbation separate in x and y . For example, eq. H21|l enforces the two eigenvalue 
equations in the bulk 



(□ - ml) Q„ = 



(24) 



Similarly, eq. H1U|I gives an eigenvalue equation for lAn ■ Combining it with the constraint equations to eliminate terms 
in l>Jj and l>" , we get 



A' 

4*' 



A^V ^ AU" 



3M^ 



2 4 y 4 

Finally, we insert the decomposition (|23|) into the action Making use of the last two equations, we get 



(25) 



5 = ^ Cnm j di^X Qra [ □ - TO^ ] Qn , 



dy Vjn Vn 



AA' 



3M^ 



yo 



dy 
A 



3AP 



(26) 



The second expression for Cmn has been obtained relating v and $ through U (cqs. (IIOII and pOfl 'l. Hermiticity and 
orthogonality of different modes are manifest from it. 

This completes the diagonalization of the starting action in terms of the different scalar modes. We have thus 
obtained the exact free actions for the physical modes 



St" = Cn j d^x Q„ [ □ - m2 ] Q„ , 



3 A/3 fvo dy 



ry° dy 3M^ / ,2 i f ^2s. y 



(27) 



The last condition determines the coupling of the scalar modes to fields on the observable brane. In the next section 
we discuss this issue in detail for a couple of specific examples. 



In the step denoted by dots we have differentiated the definition of v and made use of I12i to eliminate and 5<p' 
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IV. COUPLING OF THE PERTURBATIONS TO BRANE FIELDS 

In the previous section we have determined the canonically normahzed scalar modes Qn which enter in the exact 4 
dimensional description of the system. Their coupling to SM fields is encoded in the induced metric 7 at the observable 
brane. From the Standard Model action 

'S'sM = / d'^x^/^ CsM (28) 



we get 

5i, 



i ^ $„ (2/0) j d^xV^Qn {x) TH (x) + . . . , (29) 



where is the trace of the stress energy tensor of SM fields, while 70 denotes the background induced metric. Dots 
denote higher order couplings, as well as other possible interactions, as for example the one arising from the coupling 
of the Higgs boson to the curvature j20], which will not be considered here. Hence, the coupling A„ of the n — th 
mode is given by 



^ _ Qn rpi_i 

•'-int.n — . -L 



1 




$n (yo) 


A„ 




2 



(30) 



This shows explicitly that the normalization coefficients Cn obtained in the previous section are needed to determine 
the coupling of the modes to SM fields. 

In the following we compute the mass spectrum and the coupling of the scalar perturbations for two specific models. 
As remarked in the previous Section, we can (with a very good accuracy) neglect the impact of the brane fields on 
the wave function of the radion (they can be treated as probe fields). The first model we consider is a variation of 
the Goldberger-Wise mechanism 0], due to DeWolfe et al. (DFGK) 0, for which the background equations can be 
solved analytically. The study of the perturbations for this model was performed in details in Q in the limit of small 
backreaction of the scalar field </>. In subsection IIV Al we extend this study to large </>. The second model, discussed 
in subsection IIV Bl is a toy example which does not solve the hierarchy problem. We discuss it here because both 
the equations for the background and for the perturbations can be solved analytically. The two models show similar 
features: the mass scale of the radion is given by the backreaction of the scalar field, whereas the mass scale of the 
remaining KK modes is set by the size of the extra dimension. In addition, the coupling of the radion with SM fields 
is in both cases stronger than the one of the KK modes. 

A. The DFGK Model 

DeWolfe et al. provided a procedure to "construct" bulk and brane potentials for the scalar field for which the 
background equations can be solved analytically 7]. One specific example studied in is characterized by the 
following potentials ^ 

V = -QkHl^ +{2ku+—^ ' 



2 y ^ 6 M3 ^ ' 

U± = ±{QkM^ ~-u4)'i)T'iu4)±{(t)~(t>±)^'\^l±{(t>~(f)±f . (31) 

The upper and lower sign refer to the brane at y = and y ~ yo , respectively. The quantity 0+ (0-) denotes the 
value of the scalar field at the brane at y = {y ~ yo). The "mass parameters" fi± do not affect the background 
solutions (since they do not contribute to the background junction conditions 0). However, they enter in the 
junction conditions H12|) for the perturbations. As we remarked, the /i,± — > 00 limit forces S(f) = a,t the two branes. 
For simplicity, the following calculation is restricted to this limit. 

Exact analytical background solutions are computed in normal coordinate w for the bulk, defined as 

%(j/)dy (32) 





In the notation of 01, this corresponds to the "superpotential" W (0) = fikh/fi — ucfp 
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(the two branes are located at w = and at wq = w (t/o) ). The background solutions are 

^2 



A = exp 



-ki 



12 AP 



1) 



(33) 



where the normalization A{0) — 1 has been chosen. Following we define the parameter P = (/)^/(2M^) to 
characterize the strength of the scalar field (j). The I <C 1 limit corresponds to negligible backreaction of the scalar 
field on the background geometry. The case I = corresponds to the RS background solution JJ. 
The four dimensional Planck mass is given by 



wo roc 7i,r3 ^ 

k 



M; = 2Ivr / dw {w) ~ / dw {w) 
Jo Jo 



k k 

u ' u 



where i-Fi is a confluent hypergeometric function. For u <ti k , this expression reduces to 



(k < M) 



The electroweak scale at the second brane is generated from the warp factor 

TeV 



A (wo) 



(34) 



(35) 



(36) 



In the present model, the equations for the perturbations can be solved analytically only in the small I limit. The 
computation, summarized in Appendix^ gives at leading order in / 



Radion 
KKmodes 



1 



Arad 
1 



kyo 



2ul 



TOrad 



A„ 



V6Mp ' 
2V2IU (fcyo)'""^'' 



(kyo) 



l+u/k 



2k 



3k 



3m: k 



Mr, 



n TT 

ya 



1, 2, 



(37) 



As we mentioned, the coupling of the radion to the SM is set by the electroweak scale {A (yo) ~ ky^ ^ 10^^ — 10^^ 
sets the hierarchy between the Planck and the electroweak scales) . We also note that the radion mass is proportional to 
the backreaction parameter / , confirming that the radion corresponds to a flat direction in absence of any stabilization 
mechanism (an additional significant suppression to TOrad is given by the u {kyo) "^'^ prefactor). As a consequence, 
the radion mass is expected to be below the physical cut-off scale for the ?/o brane, set by eq. (|35|l '). On the other 
hand, the masses for the KK modes are rather insensitive to I , and are set by the inverse size (~ TeV) of the extra 
dimension. Finally, we observe that the coupling of the KK modes to brane fields vanishes in the I limit. This was 
expected, since in this limit the mode $, which is the only one directly coupled to brane fields in eq. H29|) . coincides 
with the radion. 

Eqs. (|37|) are valid for / ^ 1 . For larger /, the equations for the perturbations cannot be solved analytically. 
However, they can be solved numerically, and the coupling to SM fields can be obtained by numerical evaluation of 
the condition H27|l . which is valid for arbitrary background. We found convenient (i.e., numerics do not develop large 
hierarchies) to work in terms of the variable Y = $ , and in normal coordinate w . The bulk equation for Y is 
given in Appendix ^ It is supplemented by the junction conditions Y' = Q at the two branes. Hence, we have a 
boundary value problem, which we solve by the shooting method: in the numerical integration, we start from Y — 1 
(normalization yet to be determined) and K' = at w = , choose an arbitrary value for the mass, and evolve the 
bulk differential equation for Y . The integration gives a value of Y' {w = wq) which depends on the choice of m . 
In general, Y' (wq) ^ 0, signaling that we started from a "bad value" of the mass m. Only when Y' (wq) = we 
have found a physical mode which solves all the equations for the perturbations. This occurs only for a discrete set 
of m , which constitute the masses of the physical modes of the system. A dense scan in m allows to determine the 
spectrum of the theory, and the wave function associated to each mode. The normalization of each wave function is 
then found by numerical evaluation of (|27|l . 

Figure ^ shows the masses (left panel) and the couplings of the radion and the first (lightest) KK modes as a 
function of I . Besides I , the system is determined by the four parameter M, fc, u , and yo (or, equivalently, ). We 
fixed ujk — 1/40 (too big or too small values give too small radion masses), and k — O.lAf. The two additional 
conditions that we need are the recovery of the Plank mass in the 4D gravity, eq. (|34|l . and of the electroweak scale 
at the observable brane, eq. (|35|l . The numerical results agree with the analytical expressions (|37|l at low I. The 
radion mass grows linearly with / , while the modes of the KK modes are constant. As / increases, the masses show 
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0.1 1 10 0.1 1 10 

I I 



FIG. 1: Mass (left panel) and coupling (right panel) of the radion and the first KK modes as a function of the (rescaled) 
amplitude I of the scalar field. Parameters are chosen as fc = 0.1 Af , it = fc/4Q , and in a way to maintain the same warping for 
aU Z . 



a rather sharp decrease, which extends also to greater values of I than the ones shown in the Figure. The decrease 
is probably due to the increase of wq , necessary to preserve a strong warping as I increases (cf. eqs. ijH^ and (j^ ). 
We see that intermediate values for / are needed, since the radion mass becomes too small both at small and large I . 
The maximum is achieved for I < 10 , far from the I ^ 1 regime. Also the couplings of the modes to the observable 
sector, shown in the right panel, agree with the analytical results at small I . The coupling of the radion is constant, 
while the one of the KK modes increases linearly with I (This behavior has been explained in the paragraph after 
eq. (EH))- The coupling of all the modes sharply increases at large I , when the backreaction of (p cannot be neglected 
any longer. Both the decrease of the masses and the increase of the couplings occurring at large I increase the potential 
detectability of the bulk excitations, and stringent experimental bounds should be expected in this limit. 




0.1 1 10 0.1 1 

I I 

FIG. 2: Same as in fig.Cl but with fc = M . 



An analogous calculation is summarized in Figure where we show the masses and the couplings of the scalar 
modes for the choice of parameters k = M , u = k/AO (j/q and M are then fixed as in fig. The numerical results 
shown in the two Figures exhibit the same qualitative behavior. The increase of k is responsible for the higher values 
of the masses shown in fig.EJ The values of the couplings are instead very weakly sensitive on k , as eqs. (|37ll indicate. 
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B. A Toy Braneworld 

We consider a second example, characterized by the scalar potentials 



(see the previous subsection for the notation). 

The background solutions for this model (given in conformal coordinate y along the bulk) are 

0^M3/2[ay + /3] , A^?^%^e-^^Aoe--y/^^ , (39) 

Wo 

where a and /3 are arbitrary constants (of mass dimension 1 and , respectively). Hence, the scalar field is linear in 
y . One can arrange for a strong warping between the two branes, a S> 1 . However, it is obtained through a strong 
hierarchy (in terms of the parameters entering in ]3H\i ) between the potentials U± at the two branes. In this sense, 
the model does not constitute a natural solution to the hierarchy problem. However, it is rather instructive, since the 
bulk equations for the perturbations of this model can be solved analytically (we have z"/z = 3a^/4 in eq. H^lfl ). 

The junction conditions for the perturbations can be also solved analytically in the /i± oo limit The mass 
spectrum is 



2 /Sa'^/nTT^^ 



Wrad = \ n'^ ' "^KK = \ — h • (40) 

V 3 V 4 V / 

The mass of the radion vanishes in absence of warping, while the KK modes have a contribution proportional to the 
warping in addition to the standard KK mass ~ mr/yo (n = 1, 2, . . . ) . The masses of the radion and of the KK 
modes become comparable at large warping, ayo » 1. Finally, the coupling of the different modes to brane fields at 
the observable brane are given by 



1 



Arad 1< 



vl/2 



21/2 33/4 ^^3/2 A 



3/2 



g\/3Q 1/0/2 



1 



A 



KK 



^0 



4n7ra^/^ 1 

3M3/2yl3/2 («yo)3/2 



^V3ayo/2 ^ ^^^^ 



As in the previous example, the radion is more strongly coupled than the other KK modes, although the couplings 
have comparable strength. 



V. NUMERICS OF THE PERTURBATIONS WITH THE BRANECODE 



The analytical results derived above can be be supported by very different means, namely by the numerical study 
of the braneworld dynamics. Suppose we have a stable flat brane configuration plus fluctuations around it. This can 
be achieved by taking initial conditions which are slightly displaced from the stable configuration, so that the system 
will start oscillating around it. In this way we can study the properties and the time evolution of the perturbations, 
and see what relation they have with the eigemodes of the system computed analytically, eq. H27() . The evolution of 
the interbrane distance D (t) is often characterized by several oscillatory modes, superimposed to a lower frequency 
oscillation. Correspondingly, the Fourier transform of D (t) shows the presence of resonance bands: the one with 
lowest frequency is identified with the radion field, while the other ones with the higher KK modes. The numerical 
values of the frequencies of the bands can be compared with the analytical expression for the masses of the modes. 
In addition, we can split the time evolution in successive time intervals (each containing several oscillations), and 
perform the Fourier transform in each of them. This allows to study how the excited modes evolve in time. 



This is also obtained from the general prescription given in 01, using the "superpotential" W (0) = Wo exp U/ 1 VsmA . 

Solutions for finite large ii± , can be given as an expansion series in 1/ fi± . However, they arc not illuminating, and we will not present 

them here. 
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We make use of a numerical code, named BraneCode [l^l, designed for the study of the broad class of brane mod- 
els (O we are interested in. The numerical integration is performed with the assumption of homogeneity and isotropy 
along the brane coordinates, so that the system solved is effectively two-dimensional (time and bulk coordinates). 
This limitation removes the tensor modes from the numerical evolution (since there are no tensor modes in two dimen- 
sions), which is welcome for the present analysis. However, it also restricts the study of the scalar modes to the large 
wavelengths limit. Since we numerically integrate the exact 5D Einstein equations, the present study goes beyond 
the linear regime for the perturbations, which limits instead the analytical calculations. In the example shown below, 
the nonlinear dynamics is relevant at the earlier stages of the evolution. 

Some attention has to be paid to the choice of the initial conditions, since they have to respect two constraint 
(nondynamical) Einstein equations. The simplest possibility is to search for static bulk configurations, with 

ds^^Aiyf[-df+c^"'d^^ + dy^] , (42) 

and with (j) = (j){y) . As discussed in 0| , the class of potentials considered in section IIVI admits at most two static 
configurations. When they exist, they are characterized by different values of H , which is proportional to the Hubble 
parameter measured by brane observers. The solution with higher H is unstable, while the other one is stable. As 
explained in p^. a non- vanishing H gives a negative contribution to the physical masses of the scalar fluctuations 
around the solution, m1^^ < ttIq — 4 . In many cases m^^^. is negative, meaning that the corresponding static 
configuration is subject to a tachyonic instability [iSj. The instability typically leads towards the stable configuration, 
characterized by the lower expansion rate [T^ . 

A systematic algorithm for the search of static configurations is described in 0|. In the present discussion, we 
focus on the model H38|) . and start from an unstable static solution. Figure 13 shows the relaxation of this configuration 
towards a stable solution, of the class discussed in the previous section. We plot the time evolution of the interbrane 
distance, defined as 

ryo 

D{t)^ / dyA{t,y) , (43) 
Jo 

together with its Fourier transform. We show two different spectra, according to different (consecutive) time intervals 
immediately after the "transition" between the two solutions. 




FIG. 3: Relaxation towards a stable braneworld configuration in the model l|88^ . The evolution of the interbrane distance D (t) 
and its Fourier transform are shown. The arrows beneath the frequency axis denote the masses of the radion and of the first 
KK modes according to eq. H4()^ . Parameters are chosen as q j/o = 2, /3 = 0, e~" Aq fi+ = Ao fi- = 100 M , and the figure 
is shown in units M — 1 . 

The numerical result confirms the rapidity of the transition, and the initial excitation of the radion mode. The 
subsequent evolution can be roughly divided into two stages. During the earlier one, we can see the excitation of modes 
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of higher frequency superimposed to the radion oscillation. The positions of the modes are in excellent agreement 
with the masses obtained analytically - see eq. (|40|l - which are indicated by the arrows below the frequency axis. 
These confirm that the modes excited are precisely the ones found through the linearized analytical calculation. We 
note that the excitation of the higher modes starts already during the first oscillation of the radion field. It is also 
possible to see a decrease of the amplitude of the radion oscillations during all this first stage. The second stage is 
instead characterized by a constant pattern of oscillatory modes. The spectrum of the perturbations docs not change 
appreciably during these later times. 

These features can be readily understood. The large amplitude of the radion oscillations is due to the strong 
tachyonic instability of the first configuration. The excitation of the higher frequency modes is due to the coupling 
of the perturbations which occurs at the nonlinear level (as we mentioned, the numerical integration automatically 
takes into account the nonlinear dynamics of the system) . Most of the energy of the radion field is transferred to the 
higher modes, as also indicated by the decrease of the overall amplitude occurring in the first stage. At later times, 
the amplitudes of the oscillations become sufficiently small so that nonlinear effect can be neglected. The fact that 
the patterns of the oscillations become constant confirms that the modes are decoupled at the linear level, as eq. H27|l 
indicates. 

In the present computation, we did not include any brane fields. Their coupling to the radion may be responsible 
for some energy transfer from the bulk to fields on the brane, during the earlier stages of cosmological evolution. This 
possibility is discussed in greater detail in the discussion section. 

VI. DISCUSSION 

We considered 5D braneworld models with the inclusion of a bulk scalar field for the sake of generality and, more 
importantly, to provide stabilization of the extra dimension. Observers on a 4D brane embedded in this space will in 
general experience corrections to the 4D Einstein gravity, as well as additional four dimensional degrees of freedom 
emerging from the bulk gravity and the bulk scalar. These excitations interact with the fields of the Standard 
Model (or more generically, of the observable sector) living on our brane. Therefore, their nature and the details 
of their couplings to the observable sector arc the subject of particle physics phenomenology and experiments. The 
set-up we studied is a specific example of more general systems characterized by a higher dimensional warped bulk 
geometry, which emerges not only in the context of 5D braneworlds, but also in many other models motivated by 
string cosmology, see e.g. For this reason we will comment not only on the results of the explicit calculations of 
the paper, but also on potential extensions of our methods and results for the broad context of the warped geometries. 

In the present paper we focused on the scalar excitations of the system. The warp bulk geometry can be strongly 
curved, and the full machinery of General Relativity is required for the study of the bulk gravity/scalar dynamics. In 
general, the physical modes associated with the excitations of the coupled gravity/scalar system do neither coincide 
with the fiuctuations Scj) of the scalar field nor with the scalar perturbations $ of the metric, but rather with some 
linear combination of them. This is sharply different from the identification of the eigenmodes in the familiar KK 
dimensional reduction, where higher dimensional metric components decouple from the bulk fields. A much more 
similar context is the one of 4D cosmology with scalar fields, where cosmological scalar perturbations of the metric 
are sourced by the fiuctuations of the infiaton. In the present paper we extended this approach to the more complicated 
set-up of higher dimensional warped geometries. 

This analysis can be pursued in several directions. 

1. We computed the coupling of the radion (the lightest eigenmodc of the warped geometry excitations) and 
of the higher KK modes to Standard Model particles, which is one of the most crucial pieces of information for 
phenomenological studies. The physical modes emerge from the decomposition and diagonalization of the second order 
action for the perturbations. The diagonalization requires great care, since the starting action for the perturbations is 
particularly involved, and the identification of the relevant dynamical variables is far from being obvious. Contributions 
both from the bulk and the branes are present in this action, neither of which can be immediately diagonalized. In 
particular, a priori one could have wondered whether the two different contributions could have been simultaneously 
diagonalized, or if some mixing between the different KK levels would have unavoidably remained. 

We performed this tedious decomposition in the case of 5D braneworlds with bulk scalar field and two orbifold 
branes at the edges. The answer happens to be rather transparent and simple: the total action for the perturbations 
from the bulk and the two branes can be diagonalized in terms of free, non-interacting scalar KK excitations. The 
canonical quantization of the actions of these modes specifies the normalization of their wave function along the extra 
dimension. This in turns defines the value of the perturbations at the observable brane, which sets the coupling 
between the scalar modes and the SM fields (the overall normalization cannot be obtained from the 5D Einstein 
equations for the perturbations, which are all linear in the perturbations). For the reasons just remarked, it should 
not be a surprise that the couplings are very sensitive on the details of background bulk geometry and of the bulk 
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scalar field distribution. We illustrated this dependence with some specific examples, see for example Figs. ^ and |5J 

One of the examples we have considered is related to the Goldberger-Wise mechanism Q for the stabilization of 
the Randall-Sundrum model More specifically, we consider the model by De Wolfe et al. (DFGK) where the 
background equation for the stabilized geometry and the bulk scalar field can be solved analytically. It is convenient 
to introduce a parameter I (6(| - proportional to the amplitude of the scalar field - which characterizes the departure 
of the model from the pure Randall-Sundrum model without stabilization. Thus, I = corresponds to the Randall- 
Sundrum model, small / <C 1 to small backreaction of the background scalar on the bulk AdS geometry, while large 
I >, 1 implies a significant modification of the geometry due to the presence of . It is important to bear in mind that 
a large warping (large difference of the background conformal factor A between the two branes) can be obtained for all 
I . Thus, the virtue of the braneworld models to solve the hierarchy problem is preserved also for cases, characterized 
by large I , in which the bulk geometry may significantly differ from the AdS geometry of Q . Indeed, in the specific 
example we have considered we were able to preserve the same warping for all values of I . 

We shall make comment on the physical reliability of the models with small vs. non-small /. It is not easy to accept 
that the bulk scalar field gives only negligible, perturbative contribution to the background bulk geometry, and that 
its only role is to provide a tiny lift of the radion mode, without affecting the geometry in any other way. As we 
already remarked, the idea of RS is to avoid strong hierarchies besides the one given by the warp factor. If I is limited 
to the small values needed to remain in the perturbative regime, the radion mass turns out to be significantly smaller 
than the electroweak scale. Although this is certainly a much milder hierarchy than the one between the electroweak 
and the Plank scales, it can be avoided provided / is sufficiently large. Another point in favor of a large / is related 
to the dynamics of the radion in the early universe. To appreciate this, we computed and showed in Appendix lO 
the 4D effective potential for the radion (this effective potential is only approximately valid at small / , and should 
not be confused with the exact 4D description given in the main text). In terms of the canonically normalized field 
(corresponding to the radion in this effective description), the local minimum at which the stabilization occurs is 
minuscule. If we are interested in the early dynamics of the braneworld, it is very hard to understand how the system 
can go to such a minimum, unless it starts extremely close to it. This is a very similar problem to the one of moduli 
stabilization in hcterotic string cosmology j23| . This problem is enhanced for small I . These reasons are a strong 
motivation for studying the system also beyond the small / limit. 

The coupling of the radion to SM fields was already computed by Csaki et al. 6] in the limit of small backreaction, 
Z <C 1 . Our formula reproduces their result in this limit. In addition, we found that the strength of the interaction 
remains almost constant in the range < / <, 10. Most interestingly, however, is that for I > 10 the value of the 
coupling begins to grow sharply as a function of /. Therefore we conclude that, in the model considered, a stronger 
stabilization can enhance even by orders of magnitude the coupling of the radion to SM fields. 

Another very interesting property we found is the strong dependence of the masses of the radion and of the KK 
modes on the parameter I . In the range < / ;S 10 the mass of the radion increases linearly with / (consistently 
with the fact that the radion becomes massless as I —* 0), while the masses of the KK modes are nearly constant. 
At Z 10 the masses of all the modes are instead sharply decreasing. The decrease is probably due to the larger 
interbrane separation which must be imposed at large I - in order to preserve a strong warping at the observable 
brane - and which results in lower masses from gradients in the y direction. Hence, both very small and large / are 
phenomenologically excluded, since they lead to a too light radion. The radion mass is maximized for I ^ 10 , far from 
the perturbative / <^ 1 regime. This paper provides the technical tools for extending the phenomenological studies to 
this relevant regime. 

2. We identified the KK eigenmodes of the braneworld excitations in the general case, without any assumption of 
small backreaction of the bulk scalar field. In the linear regime they are free non-interacting dynamical fields, since 
they correspond to canonically normalized eigenmodes of the action expanded up to second order in the perturbations. 
However, we may expect that they are coupled at higher (third and further) order in the decomposition of the full 
action. To verify this, we made use of numerical simulations of the dynamical evolution of the braneworlds, which 
give complementary informations to the analytical study. The simulations are based on the BraneCode, which was 
developed to solve the fully nonlinear self-consistent evolution of the braneworld with the scalar field . We started 
with an unstable brane configuration, which evolves to a stable one plus excitations around it. The excitations are 
composed of a superposition of oscillatory modes, whose frequencies (determined by a Fourier decomposition) coincide 
with the eigenmasses found by the analytical study of the linear regime. The higher KK modes are excited by the 
nonlinear interactions with the radion field during the earlier stages of the evolution. At later times, the nonlinear 
interactions can be neglected, and the pattern of oscillations becomes constant. 

The discussion of the brane system restructuring and of the excitation of its eigenmodes leads us to another subject 
where the interaction of the radion and the other KK modes with SM particles is important. This is (p)reheating of 
the universe after inflation, and, in particular, the role that the brane excitations may have had in the generation of 
brane fields. Inflation in the braneworld setting corresponds to the curved brane, so that after inflation the brane is 
flattening as the brane curvature decreases with the decrease of the Hubble parameter H{t). As we just noted, the 
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reconfiguration of the brane geometry leads to interbrane oscillations, which correspond to excitations of the radion 
and the other KK modes. The interaction of the bulk modes with the SM particles will results in the decay of these 
oscillations. More specifically, let us consider a coherently oscillating radion (f> = $o cos nirad t interacting with some 
scalar field x of mass m living on our brane. The equation for the amplitudes of the eigenmodes Xk{t)e^^^ of quantum 
fluctuations of x can be reduced to 

Xfe + ^fc^ + - (m^ + ^TO™d)*o cos m™d Xk ^ (44) 

(more accurately, Xk refers to the canonically normalized field, and m to the physical mass, both obtained after a 
rescaling with the conformal factor at the brane; in addition, one should take into account that the amplitude <i>o of 
the oscillations is decreasing due to the cosmological expansion of the brane). This is the familiar equation for the 
parametric resonant amplification of quantum fluctuations interacting with the oscillating background field (25^ . The 

strength of the effect is given by the parameter q = (-^1 h ■^) $o- The amplitude of the radion oscillations (here 

$0) corresponds to the induced metric fluctuation, which is below unity. Thus, in case of a mass ratio — < 1 , 

2 

the parametric resonance is not strong enough for preheating. In the opposite case, — 3> 1 , the parameter q is 

9 2 

large, but the effect works only in the higher instability bands where k > — , and it is again very small. We 

conclude that it is hard to have preheating effects from the radion oscillations, so that the radion oscillations decay 
in the perturbative regime. An interesting earlier paper [26| seems to conclude that preheating effects associated to 
the oscillations of the radion are instead more significant. We think that, once the expansion of the brane. and the 
consequent decrease of the amplitude of the oscillations, are taken into account, the result claimed in |2g may be 
significantly reduced. 

3. Finally, we will discuss other virtues of the method of the full action series decomposition which we used in the 
paper. While we focus on the problem of identification of the eigenmodes of the gravity/scalar dynamical system in 
order to quantize the radion and the higher KK modes, in principle we can extend this method to discuss another 
important sector of the effective 4D theory, namely effective 4D gravity. There is extensive literature on the properties 
of the 4D gravity in the braneworld models; here we discuss only the low energy limit, that is the linearized gravity at 
the brane |17ll2a . l27l | . On general grounds, one expects that the 4D effective gravity at the brane will be a Brans-Dicke 
(BD) theory. For pure RS model without brane stabilizations, the BD gravity contains a massless scalar, which is the 
radion. For stabilized branes the radion acquires a mass, so that one expects a BD theory with a massive BD scalar. 
The latter situation was studied in 28], where the Lagrangian for the 4D gravity was derived as a series containing 
the Einstein term R plus higher derivatives term of the form /tti^^^, plus even higher derivatives terms. 

The main method adopted in the literature for the derivation of the low energy 4D effective gravity is based on 
the introduction of a point-like stress-energy source at the visible brane, and on the investigation of the gravitational 
response upon this source, including the bulk gravity, the scalar field, and the shift in the embedding of the branes \2T\. 
Notice, however, that the signature of the low energy gravity (BD or R^ theory) is manifested in the free gravitational 
modes even without the stress-energy sources. Therefore, in principle we can address the question regarding brane 
gravity in terms of the series decomposition of the full action. To do so, we need to include not only the scalar metric 
perturbations $ in the form lO, but also the transverse-traceless gravitational modes /i^iB- Bulk gravitational waves 
in general have 5 degrees of freedom. However, the massless scalar and vector projections of the bulk graviton are 
absent, so that the massless graviton has only 2 degrees which corresponds to usual 3-1-1 dimensional transversal and 
traceless (TT) gravitational waves ft^J . We can extend 10 to the form where the full 5D metric perturbations 
Ham are 

/i55 = 2$ , /i5^ = , h^^ = /i^J - 77^^$ . (45) 

We then have to include TT modes in the decomposition of section Fill Al Up to the second order in perturbation 
series, the effective 4D action acquires the new term hJ^'^.phJ^J'-P /A which is equal to the linearized 4D curvature R, 
see e.g. 0|. The BD coupling (1 -|- Q)R emerges only in the third order of the perturbation series, which is beyond 
the scope of this paper. Here we just note that a terms with the structure $^JJ'''^^Jp will appear in the third order 
decomposition. 

What is important, is that with the present approach we explicitly see that the BD scalar is a massive one. This 
follows from the second order decomposition of section UlI Bl Now, how the BD theory with the massive radion (as a 
BD scalar) can be reconciled with the R + R^ theory derived in |28||? It is well known |29(| that a theory of gravity 
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with the higher derivatives in the form £1 = ^ igf yR + gm^ j conformally equivalent to the Einstein theory with 



the scalar field (f) and with Lagrangian £2 — ~t^R + h4'tJ.4'^ ~ ^(0)- The metrics of theories Ci and £2 are related 
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by conformal transformation g^^ (l 



R 

3m ^ 



g^^. The scalar fields 4> is obtained from the scalar curvature R as 
R 



3M2 

^ In 1 - 

167r 



and the potential is given by formula 



v{4>) 



327r 



exp 




(46) 



(47) 



In the limit of small cf) we obtain V{(f)) w ^rn^cj)^. Thus, for small values of the radion (i.e. BD scalar) the Einstein 
gravity with small high derivative corrections is equivalent to the theory without R^ term but with instead the 
additional massive scalar. 
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APPENDIX A: COMPUTATION OF S2 

In this Appendix we outline the computation of the action (|18|l for the perturbations. Following |19| |. we first 
perform an ADM decomposition along y of the gravitational part of the action One can easily verify that, for a 
line element of the form 



£■ {t, X, y) {-dt^ + dy?) + {t, X, y) dy^ 



(Al) 



one has 



\G\R^ y/\G\iR^ 



£(4) 



12 A'^ 





'A^ A'' 







2f]^''d,,{A^d,B) 



(A2) 



where i?(4) ~ —677'"' {df^di^A) j A? is the curvature of the four dimensional slices. The last term in (|A2(I is a boundary 
term on the 3 + 1 slices, and it can be dropped. On the other hand, the total derivative in y does not vanish, since 
the space in that direction is limited by the two branes. However, its contribution to the action precisely cancels with 
the Gibbons-Hawking term -y— 7 \K\ in the starting action © for the two branes. 

Hence, only the first term of (|A2p is relevant. Combining it with the bulk action for the scalar field, we get 



•Sbull 



Va 



d^xA^B 



-6 T] 



3 „tii^]2L 



dud.A UAPA'^ 



A^ 



S2 A'^ 



2A^ 
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2 62 



12 
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(A3) 



The coefficients A and B are obtained from the line element ^ . Expanding the bulk action at second order in the 
perturbations, and making use of the background equations to simplify the final expression, we obtain the bulk part 
of the action H18(l . The brane contributions to 118|l are instead the expansion of the last term of 



APPENDIX B: PERTURBATIONS IN THE DFGK MODEL FOR I < 1 



In the following, we derive the analytical results given in H37|l . valid in the I <C 1 regime. Some of the present results 
were also obtained in . It is useful to work in terms of the rescaled mode Yn = A"^ , and in normal coordinate 
w . By combining the equations for the background quantities and the perturbations, one finds 

Yn = , f,:U=o,»o=0 (Bl) 

(in this Appendix, prime denotes differentiation with respect to normal coordinate w ) . 



y:'-2{^ + ^]y: 



2^-2. 
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For the moment, we restrict the computation to Z = , that is we ignore the backreaction of the scalar field on the 
background geometry (notice that (j)" /(j)' = —u is regular in this limit). The bulk equation is then solved by 



with arbitrary constants y}P^ and C„ . Differentiating this solution gives 



y(0) g-«.(fc + u) 



Y' 



-fiir, e- 



J. 



2+u/k 



(^e'=-)+C„y2+uA. 



The boundary condition at w = can be used to enforce 

J2+u/k (t-) , 



c 



Y2+u,k{^) r(2 + f)r(3 + 



\2k) 



(B2) 



(B3) 



(B4) 



where the last limit holds for rrin/k — * . As we will now see, this is the case we are interested in, since m is set to 
the electroweak scale, while k ^ Mp . 



Eq. I|B3|) shows that a mode ljB2|l with vanishing mass solves also the boundary conditions. This mode corresponds 
to the radion, whose mass is due to the backreaction of on the geometry: a non-vanishing radion mass is found 
only once the bulk equation (|Bip is expanded at the first nontrivial order in I . Let us postpone this calculation till 
the end of this Appendix, and concentrate on the KK modes, for which the present calculation (/ = 0) is accurate 
enough. As clear from the expansion (jB4|) . the second term in ljB2p is completely irrelevant for m„/fc ^ TeV/Mp. 
The junction condition at w = wo then gives 



J (rUn 



k Wo 







(B5) 



Hence, the masses of the KK modes are determined by the poles of this function. By expanding it for large argument 
we find 



TT 

2/0 



15 + 16m/A; + 4MVfc2 



2k 



n= 1, 2, 3, 4, 5,. 



(B6) 



This approximation becomes more and more accurate as n increases. However, the comparison with numerical results 
shows very good agreement at all n . The normalization coefficient Y^'^ is computed by inserting the solution l|B2p 
into the normalization condition (|27|l . Since we are interested in the leading result for small Z, only the first term 
in lf?7|) has to be computed. We get 



. ^2kw j2 fm-n kw\ 

dwe ^2+t(— e j 



16 k P V? 



2 „2kwo 
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-2 fmn 
(— ' 



k wq 



(B7) 



We can now evaluate the normalized $„ = Yn/A^ at the brane. We obtain, at leading order, the coupling of the 
KK modes with brane fields reported in eq. H37I) . 

We still have to compute the normalization and the mass of the radion field. As remarked, the mass can be obtained 
by inserting the background solutions H33|l into the bulk eq. ljBl|l , and expanding at the first nontrivial order in I . 
The normalization can however be found already from the / = term. At zero-th order in I , the radion is massless, 
and eqs. (|B1|I are solved by a constant Y (as it is also confirmed by the limit of (|B3|I and l|B4p for rn„ ^0). We 
thus recover the result $ cx A~'^ given in H14|) . Inserting a constant Y in H27|) we find the radion normalization, and 
finally the coupling to SM fields given in eq. H37|l . The radion mass is obtained by expansion of the bulk eq. ljBl|) at 
second order in I , using the ansatz (Yq is the constant 0— th order solution) 

Y {w) = Fo + l^y [w) , = Pfri^ . (B8) 

This equation can be easily solved. Imposing that y' (w) vanishes at the two branes gives the radion mass 

^ 4/2^2 (2fc + u) (e2 1) 4Z2u2(2A: + u)e-2('=+")«'o 



m 



3 fc (e2(2fc+")"'o _ 1) 



3fc 



(B9) 



The function y (w) modifies the radion normalization at subleading order in / ; the effect can be neglected provided I 
is sufficiently small. 
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APPENDIX C: EFFECTIVE POTENTIAL FOR THE RADION 



In this Appendix we compute the effective 4D Lagrangian for the radion in the model (|31|1 . following the original 
computation of The 4D Lagrangian is obtained starting from an ansatz which does not satisfy the 5D Einstein 
equations of the theory, so it should be considered only as a first, heuristic study of the system. We discuss it here 
mainly to compare it with the exact 4D description obtained in section IIIII To derive it, it is convenient to use 
"rescaled" normal coordinates 



A (w)^ Tj^jy dx^ dx'^ + wq {ty^ dw^ 



w G [0, 1] 



(CI) 



so that the interbrane distance is encoded in the (time-dependent) parameter wq ■ The advantage of this parameteri- 
zation is that the kinetic term for wq is very easily obtained. As in we neglect the backreaction of the scalar field 
(j) on the metric; hence, the analysis is restricted to the small / limit (in addition, the effective Lagrangian only refers 
to the radion field, while the result (|27|l is valid for all the physical modes). We "promote" the parameter wq to be 
time-dependent, so that 



— g-fcu)u)o(t) 

The kinetic term for wq is obtained by integration over w of the i?(4) term in ljA2|l . which gives 

I 2 
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dw wqA 
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A3 
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-kwa{t) 
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(C3) 



where 7 is the canonically normalized field (a further contribution to Ck is given by the kinetic term of (f> , but it is 
subdominant in / ) . 

The effective potential for the radion is obtained by integrating the action over w . The 0(^°) part cancels, 
corresponding to the cancellation of the 4D cosmological constant in the RS model. Hence, we effectively integrate 
the action for the scalar field, after subtracting the bulk cosmological constant and the brane tensions from the 
bulk/brane potentials 



cff 







/ dw A'^wq 









A{0fUa-A{\fu^ 



(C4) 



(in this Appendix, prime denotes differentiation with respect to w). The rescaled potentials are (cf. eqs. (I31II ) 

,2 " 



V = [2ku+ y 



2V2IU A/3/2 (0) _ V2 / M'^l 



2 



(0(O)-^/2;Af3/2) , 

^10 +2V2luM^/^e-""^° V2?M3/2e-"«") + ^ (1) - y2ZM3/2 e"""")^ . (C5) 



Contrary to this Appendix, in the other parts of this paper the background solution is fixed to the static configuration 
determined by the bulk/brane potentials. In the present calculation, this would amount in setting wq = qo (this can 
be easily seen by comparing the value of A at the second brane with the potential on that brane, as given here and in 
section llV All . However, our aim here is to compute the effective potential for arbitrary wo , so that the bulk ansatz 
for A (w) and for {w) has to be given also for wq qo . As we shall see, the effective potential Vcs will be precisely 
minimized for wq = qo , in agreement with the exact equations of the 5D theory. 
The bulk equation for in the unperturbed background HC2|I is 



0" - 4 fc Wo 0' - Wo (4 fc M -I- u^) = . 
In addition, we have the boundary conditions at the two branes, which, in the large /z± limit, read 

(/)(0) = \/2;m3/2 ^ 0(1) = \/2/Af3/2g-«9o , 

These equations give 
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We note that, for wq = qo , the expression IjCSp reduces to (j) = \f2l M'^/^ exp {—uq^w) , in agreement with eq. H33|l . 
The effective potential for the radion is now easily computed to be 



4 ;2 m3 {2k + u) [e"("'°-'fi) - l] 



g2(2fc+ti)tuo _ 2_ 

which is indeed minimized (and vanishes) for wq = qq . 



(C9) 
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FIG. 4: Effective 4D potential as a function of the canonically normalized field 7 . We fixed k — M, u — M/40 . The parameter 
qo is chosen so that exp (—A: go) ~ 10~^® sets the electroweak scale at the observable brane. Note that the interbrane distance 
decreases as 7 increases. 



To summarize the computation, the effective 4D Lagrangian for the radion is, in terms of the canonically normalized 
field 7 and at leading order in I , 



V^s^Al^AP {2k + u) 



-u/k 



2{2k-\-u)/k 



(CIO) 



In Figure 0] we show Vdf (7) for a specific choice of the parameters compatible with the solution of the hierarchy 
problem (see the main text). We note that the interbrane distance decreases as 7 increases. The limit 7 — > 
{wq — > 00) corresponds to infinite distance between the two branes, while the two branes coincide for 7 — ^JOM^/k 
{wq = 0). The potential Vcn sharply increases as the two branes approach each other, due to the strong gradient 
energy of <j) . On the contrary, vanishes at very large distances (small 7 ) . In addition, Vcs has a minimum at 
finite 7 = 7o , corresponding to = 90 • This is the value of the radion set by the stabilization mechanism. 

We should caution against the use of this potential for a strongly time-dependent 7 {t) , and for 7 too far from the 
minimum 70 . Indeed, the starting ansatz (|C2p and (|C8p does not satisfy the Einstein equations of the 5D theory even 



in the small I limit. More accurate answers can be obtained from the (numerical) integration of the equations of the 
exact 5D theory Still, the effective potential can be used to gain a qualitative understanding of the system. We 
see that 7 has to be always very close to 70 , if we want the stabilization to take place. We note the close analogy 
with to the problem of stabilizing the dilaton of superstring theory from gaugino condensation [22| : starting from a 
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static configuration far from 70 , the two branes will move to infinite distance, since the barrier in the potential at 
7 ;S 7o is too small to "stop" the evolution of the two branes, unless 7 ~ 70 from the very beginning. Although this 
problem is present at all I , the scale of the potential shows that it becomes worse as I decreases, since one can start 
from smaller interbrane distance without Vcs to exceed the limit of validity of the computation (e.g. VcS smaller than 
the scale set by the AdS curvature). 

The effective potential (|(]10(l allows a very accurate determination of the mass of the radion in the small I limit 

2 _d^Vcs, _ Afu^ {2k + u) e^'"J° 4 ^2 + e-2(fe+")'?o 
"^'^ = '7=70 - 3 fc (e2(2fc+")9o - 1) " 3fc ' ^ ' 

where the last expression holds for strong warping. We note the remarkable agreement with the 5D calculation, 
eq. (IB9|) . This result contradicts some claims in the literature of a discrepancy between the two cornputations for the 
radion mass in the DFGK model. The existence of different answers was reported for example in [6j, end of section 
6. The origin of the discrepancy is most probably due to a miscalculation of Kff in the analysis quoted by j^: the 
mass ()C11|I obtained from the effective potential computed here is in excellent agreement with the one from the 5D 
equations given in IjBQp and in ||6|. 
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